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Abstract 

In engineering applications, one of the major challenges today is to develop reliable and 
robust control algorithms for complex networked systems. Controllability and observability of 
such systems play a crucial role in the design process. The underlying network structure may 
contain symmetries - caused for example by the coupling of identical building blocks - and 
these symmetries lead to repeated eigenvalues in a generic way. This complicates the design 
of controllers since repeated eigenvalues might decrease the controllability of the system. In 
this paper, we will analyze the relationship between the controllability and observability of 
complex networked systems and graph symmetries using results from representation theory. 
Furthermore, we will propose an algorithm to compute sparse input and output matrices based 
on projections onto the isotypic components. We will illustrate our results with the aid of two 
guiding examples, a network with D4 symmetry and the Petersen graph. 

This paper is dedicated to Marty Golubitsky on the occasion of his 70**^ birthday. 

1 Introduction 

In several real world applications, the crucial task is to design reliable control and sensing tech¬ 
niques for complex networked dynamical systems. Examples of such networks are, to name but 
a few, electric circuits, power grids, buildings, multi-agent systems, consensus problems, or social 
networks. These complex systems may exhibit symmetries in a variety of different ways due to the 
construction of the network or due to intrinsic characteristic properties of the system itself. It is well 
known that equivariance properties lead to the existence of multiple eigenvalues in a generic way. 
In bifurcation theory this has first been recognized by Marty Golubitsky and Ian Stewart (cf. their 
ground-breaking book [1]). But also in control theory this fact has already been realized, see for 
example [2]. Also in the recently published paper [3], it is shown that symmetries in a network 
might decrease controllability and observability depending on the specific symmetry type. In [4], 
the authors analyze the relationship between the controllability of a networked system and graph 
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symmetries by utilizing so-called signed fractional automorphisms, a relaxation of the definition of 
an automorphism allowing not only permutation matrices but also matrices whose row and column 
sums are equal to one. The controllability and observability of networked linear systems has also 
been analyzed in [5]. However, the related characterization is based on the extension of the notion 
of strict system equivalence to networks, and equivariance properties are not directly taken into 
account. 

In this paper, we will exploit results from representation theory to determine whether a given 
system is controllable and observable and derive necessary and sufficient conditions. The idea of 
representation theory is to formulate abstract algebra problems as linear algebra problems. Each 
element of a group, for instance, can then be represented by an invertible matrix. Using this 
approach, it is possible to compute projections onto invariant subspaces and to decompose the 
problem into smaller subproblems, resulting in a block diagonal structure. Additionally, we will 
propose a new algorithm to compute sparse input and output matrices B and C for a given control 
problem. This is an important problem, in particular when only a limited number of control and 
observation input nodes is available [3]. 

The paper is organized as follows: In Section 2, the basic concepts of equivariance and represen¬ 
tation theory will be introduced. Section 3 shows how results from representation theory can be 
used to determine controllability and observability of networked systems with inherent symmetries. 
The results will be illustrated with the aid of guiding examples. Section 4 concludes with a brief 
summary and possible future work. 

2 Basic Concepts of Equivariance and Representations 

In this section, we will review the basic concepts of equivariant dynamical systems theory which will 
be relevant for the considerations within this article. For background material and a description of 
the mathematical details, the reader is referred to the classical textbooks [1, 6]. 

For simplicity, we initially restrict our attention to (uncontrolled) linear dynamical systems of 
the form 

X = Ax, (1) 

where x G M"" and A G M”’”'. Controllability and observability will be analyzed in the following 
section. 

2.1 Equivariance 

In the application we have in mind, the network structure of the underlying system induces certain 
structural properties of the matrix A. In mathematical terms, this can be expressed by the following 
equivariance condition: We assume that 


yA = Ay 

for all y G r, where T C 0(n) is a finite group of orthogonal matrices. 


( 2 ) 
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Example 2.1. The matrix 


A = 


(B Cl 0 C 2 \ 

C 2 B Cl 0 
0 C 2 B Cl 
\Ci 0 C 2 B / 


with B, Ci,C 2 G H 
the group element 


is equivariant with respect to the action of the cyclic group Z 4 generated by 


Ri = 


/O 

0 

0 


Id 

0 

0 

0 


0 

Id 

0 

0 


o\ 

0 

Id 

0 / 


where here and in what follows Id denotes the d-dimensional identity matrix. 

In fact, the matrix A reflects the structural properties of the network illustrated in Figure 1(a). 
This network is invariant under rotations by 90, 180 and 270 degrees. Note that if the coupling 
structure to the nearest neighbors is identical, then the network is also invariant under reflections 
as shown in Figure 1(b). In that case, Ci = C 2 and the matrix A is equivariant with respect to the 
group D 4 generated by the group elements 


Ri = 


In general, the dihedral group Dk consists of 2 k elements: k rotations Ro{= In), Ri, ■ ■ ■, Rk-i and 
k reflections Si,... ,Sk. A 


/O 

Id 

0 

0\ 



0 

Id 


0 

0 

Id 

0 


0 

Id 

0 

0 

0 

0 

0 

Id 

CLTld — 

Id 

0 

0 

0 

Vd 

0 

0 

0/ 


Vo 

0 

0 

IdJ 




Figure 1: (a) Network representing the coupling structure corresponding to the equivariance prop¬ 
erty of the matrix A defined in Example 2.1. The different colors of the arrows indicate 
the difference in the strength of the coupling to the nearest neighbors. These correspond 
to the coupling blocks Ci and C 2 in the matrix A. (b) 114 -symmetric network. The two 
axes of reflection are illustrated by dashed lines. 
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Example 2.2. Consider the Petersen graph V shown 

in Figure 2, 

cf. [7]. The system matrix has 

the following structure 














c 

0 

0 

C 

C 

0 

0 

0 




C 

B 

c 

0 

0 

0 

c 

0 

0 

0 



0 

C 

B 

c 

0 

0 

0 

c 

0 

0 



0 

0 

C 

B 

C 

0 

0 

0 

c 

0 


A _ 

c 

0 

0 

c 

B 

0 

0 

0 

0 

c 


A. — 

c 

0 

0 

0 

0 

B 

0 

c 

c 

0 



0 

c 

0 

0 

0 

0 

B 

0 

c 

c 



0 

0 

c 

0 

0 

C 

0 

B 

0 

c 



0 

0 

0 

c 

0 

C 

c 

0 

B 

0 



\0 

0 

0 

0 

c 

0 

c 

c 

0 

Bj 


Note that each vertex could again represent a 

d-dimensional dynamical system. However, for sim- 

plicity, we will assume here that d 

= 

1. 

The automorphism group of the Petersen graph is iso- 


morphic to S 5 , i.e. Aut('P) = S 5 . That is, there are 5! = 120 automorphisms. The group S 5 is 
generated by the transposition irt = (1 2) and the n-cycle tTc = (1 2 3 4 5). Every permutation 

TT £ S 5 acts as follows on the graph V: Each vertex with the assigned set {a,b} is mapped to the 

vertex with set { 7 r(a),' 7 r( 6 )}. Thus, induces (3 7)(4 10)(8 9), a transformation of the graph that 
fixes vertices 1,2, 5 ,6 and exchanges 3 •(-)• 7, 4 -H- 10, and 8 -H- 9. The permutation -Kc induces 

(1 4 2 5 3) (6 9 7 10 8 ), a rotation of the graph by 144 degrees. A 

2.2 Some Representation Theory 

We now briefly review some elementary concepts from group representation theory that are relevant 
within our framework. As above, let F C 0(n) be a finite group of orthogonal matrices. 

Definition 2.3. A suhspace V C M"" is called F-invariant if jv £ V for all v £ V and all 7 G F. 
A T-invariant subspace V is F-irreducible if it does not contain a proper nontrivial T-invariant 
subspace. 

Observe that the operation of F on a F-invariant subspace V induces a corresponding represen¬ 
tation 

: F ^ L(y), 7 71 ^ (3) 

of the group on that space. Here L(H) denotes the space of linear maps from V into V. Repre¬ 
sentations of F induced by irreducible subspaces are called irreducible representations of F. Up to 
isomorphism, a finite group possesses finitely many different irreducible representations. In fact, 
the irreducible representations of the finite groups are well known. 

The character xu of an irreducible representation d assigns to each element of d its trace, that is 

Xuil) = trace(r 9 ( 7 )) (4) 


for all 7 G F. 
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Figure 2 : Petersen graph V. The vertices of the graph are labeled 1,2,..., 10. In addition to 
the vertex number, each vertex is assigned one of the subsets {a,b} C {1,2,3,4,5}. 
Vertex i is then connected to vertex j if and only if the corresponding sets are disjoint, 
i.e. {ai,bi} n {aj,bj} = 0. 

Example 2 . 4 . In our guiding example - the dihedral group P = D 4 acting on M”, where n = 
8 (i.e. set d = 2 in Example 2.1) - there exist three different types of irreducible subspaces or 
representations, respectively. Here and in what follows ... ,w^) denotes the span of the 

vectors ... ,w^ G 

(a) ill : r — >• is given by = 1 for all 7 G P. ^4 corresponding irreducible subspace 

Vi C M” is given by (v^) with = (1,0,1,0,1, 0,1, 0)'^. 

(h) : T —)• ig given tjy 112 ( 7 ) = 1 for all 7 G {Rq, R 2 , Si, S 3 } and 112 ( 7 ) = —1 for all 

7 G {i?i, i?3, 52 , 5 ' 4 }. A corresponding irreducible subspace V 2 C ME is given by (u^) with 

= ( 1 , 0 ,- 1 , 0 , 1 , 0 ,- 1 , 0 )'^. 

(c) ila : r —)• ig generated by 

" 0 ^) and MSi)={j^ 

A corresponding irreducible subspace V 3 C M” is given by (u^,?;^) with 

uf = (0,0,1,0, 0,0,-1,0)^ and uf = (-1, 0, 0, 0,1,0,0,0)^. A 
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Remark 2.5. (a) We remark that the group D 4 possesses two additional irreducible representa¬ 
tion, namely '!?45 : F —>■ given by 


= 1 for all 7 G {Rq, Ri, R2, R3}, 'dii'j) = -1 for all 7 G {5i, S2, S3, S' 4 } 


and 


1 ^ 5 ( 7 ) = 1 for all 7 G {-Ro, R 2 , S 2 , 5 ' 4 }, idbil) = “1 for all 7 G {-^ 1 ,^ 3 , -S*!, S 3 }. 


In our network example, these are not relevant since corresponding irreducible subspaces do 
not exist in M"' in this case. To see this, consider the representation Here a vector 
u G M” \ {0} satisfying 

-yv = v for all 7 G {Rq, Ri, R 2 , R 3 } 

has already full D^-symmetry and therefore cannot satisfy '-yv = —v for all 7 G {5i, S 2 ,S 3 , S' 4 }. 
The analogue argument is valid for the irreducible representation 

(b) For the cyclic case Ci ^ C 2 - that is, F = Z 4 acting on M” as shown in Figure 1(a) -, the 
four one-dimensional complex irreducible representations are generated by 


MRi) = 1 , 

Over the real numbers, the 
cyclic group Z 4 is 


'd 2 {R\) = — 1 ) '&3{Ri) = F and i? 4 (i?i) = —i. 

corresponding two-dimensional irreducible representation of the 


MRi) = 


0 -1 
1 0 


Observe that eigenspaces of the matrix A are F-invariant. This follows immediately from the 
equivariance of A. That is, let v be an eigenvector of the matrix A with corresponding eigenvalue 
A, then 

Av = Xv =i> A{'yv) = 'yAv = X{'yv). 

Thus, the existence of higher-dimensional irreducible representations naturally leads to the existence 
of repeated eigenvalues of the matrix A. It is this observation which is of particular relevance for 
control and observation purposes and we will now further elaborate on this. 


Definition 2.6. (a) Let il be an irreducible representation of T of dimension n^. Then D is 

absolutely irreducible if the only linear maps on M”’’ which commute with all 'j € T are 
multiples of the identity. 


(b) A representation Dw ofT is called F-simple if either 
(i) Hw is irreducible but not absolutely irreducible, or 

(a) ilw is the sum of two isomorphic absolutely irreducible representations dy. 

Example 2.7. The irreducible representation 1)3 of D 4 is absolutely irreducible. The irreducible 
representation d 3 of Z 4 is not absolutely irreducible. To see this, note that Z 4 is a commutative 
group. A 
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Using this terminology, we can now describe the generic structure of eigenspaces of equivariant 
matrices. This result is due to M. Golubitsky and I. Stewart and can be found in [1]. 

Proposition 2.8. Let A be an eigenvalue of a matrix A satisfying the equivariance condition (2). 

(i) Suppose that A G M and let W be the corresponding eigenspace. Then, generically, W is 
absolutely irreducible. 

(a) Suppose that A G C \ M and let W C be the space spanned by real and imaginary parts of 
corresponding eigenvectors. Then, generically, W is T-simple. 

Observe that the first part of Proposition 2.8 implies that, generically, there will be no multiple 
eigenvalues (real or complex) in problems which possess Z^-symmetry. This follows from the fact 
that the two-dimensional representation of Z^ is not absolutely irreducible. However, we do expect 
to obtain multiple eigenvalues (real or complex) in networks which are D^-symmetric or which have 
the structure of the Petersen graph. 

Definition 2.9. Let d be an irreducible representation. Then we denote by its isotypic compo¬ 
nent, that is, the subspace o/M” which consists of the sum of all irreducible subspaces isomorphic 
to d. 


We are now in the position to state a fundamental result from group representation theory which 
is of importance for our investigations. 

Theorem 2.10. We have 

= (5) 

where the sum is taken over all irreducible representations of P acting on M”. Moreover, each 
isotypic component is an invariant subspace for the matrix A: 


A{W^)) C W,,. 


( 6 ) 


In representation theory, the decomposition (5) is called the isotypic decomposition. Conse¬ 
quently, A possesses a block diagonal structure if we choose a basis according to this decomposition 
of M”. 

The projections onto the isotypic components are explicitly known. In fact, let d be an irreducible 
representation of dimension n.Q. Then the projection onto the corresponding isotypic component is 
given by (see (4)) 

^ ^ Xvilhx. (7) 

' ' 7er 


Example 2.11. First we compute the projections Pi onto the isotypic components Wi for the 
irreducible representations di, i = 1,2,3, of our guiding example. We obtain by using (7) 



(Id 

Id 

Id 

Id\ 

1 

Id 

Id 

Id 

Id 

4 

Id 

Id 

Id 

Id 


Vd 

Id 

Id 

Id) 



/ Id 

-Id 

Id 

-Id\ 

1 

-Id 

Id 

-Id 

Id 

4 

Id 

-Id 

Id 

-Id 


\-Id 

Id 

-Id 

Id) 


and 
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Id 0 -Id o\ 

0 /d 0 -Id 

-Id 0 Id O’ 

0 -Id 0 Id/ 

The corresponding isotypic components are 

Wi = { (a, a, a, a)'^ € MT : a G } , 

W2 = { {a, -a, a, -aG M” : a G } , 

VF3 = {(a, b,-a,-b)^GW^:a,bGR‘^} 



For example, for n = 8, d = 2 , and 


/lO 

-10 

6 

3 

0 

0 

6 


3 

30 

1 

5 

0 

0 

1 

5 

6 

3 

10 

-10 

6 

3 

0 

0 

1 

5 

3 

30 

1 

5 

0 

0 

0 

0 

6 

3 

10 

-10 

6 

3 

0 

0 

1 

5 

3 

30 

1 

5 

6 

3 

0 

0 

6 

3 

10 

-10 

Vl 

5 

0 

0 

1 

5 

3 

30 / 


we obtain with respect to a basis, composed from bases ofWi, W2, and W3, 


/ 0 \ / 1 \ / 0 \ / 1 \ /- 2 \ / / 0 \ 


0 

1 

0 

1 

0 

\0j 


1 

0 

1 

0 

1 

0 

Vi/ 


0 

-1 

0 

1 

0 

\0j 


1 

0 

-1 

0 

1 

0 

V-1/ 


0 

1 

0 

-1 

0 


Wi W2 

the following block diagonal structure for A 


1 

0 

1 

2 

-1 

0 

V-1/ 


0 

0 

0 

-1 

0 

0 

V 0/ 


W3 


0 

1 

-1 

0 

0 

-1 

V 1/ 


/22 

-4 

0 

0 

0 

0 

0 


5 

40 

0 

0 

0 

0 

0 

0 

0 

0 

-2 

-16 

0 

0 

0 

0 

0 

0 

1 

20 

0 

0 

0 

0 

0 

0 

0 

0 

10 

-4 

-3 

-7 

0 

0 

0 

0 

3 

24 

3 

0 

0 

0 

0 

0 

6 

22 

19 

7 

Vo 

0 

0 

0 

0 

-6 

3 

27 / 
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Remark 2.12. If Ci ^ C 2 in our guiding example, we obtain the complex isotypic components 


Wi = { (a, a, a, a)^ G C"" : a G } , 

—n, a, —a)^ G C"" : a G } , 

= { («) —ia, —a, iay G C” : a G } , 

Wi = { (a, ia, —a, —ia)'^ G C” : a G } . 

Example 2 . 13 . For the Petersen graph introduced in Example 2.2, we obtain the irreducible rep¬ 
resentations: 


(a) : r —)• is again given by i9i(7) = 1 for all 7 G T. 


(b) : r —)> is generated by 




0 

0 

1 

0 \ 


/ 

0 

-1 

-1 

0 

- i \ 


0 

0 

0 

0 

1 



-1 

1 

1 

-1 

1 


0 

0 

1 

0 

0 

and i? 2 ( 7 !'c) = 


1 

0 

-1 

0 

-1 


1 

0 

0 

0 

0 



0 

1 

0 

0 

0 


\0 

1 

0 

0 

oj 


V 

0 

-1 

-1 

1 

0/ 


(c) : r —)• is generated by 




/-I 

0 

0 

o \ 

-1 

1 

0 

0 

0 

0 

1 

0 

V-I 

0 

0 

1/ 


and 'i ?3 (vTc) 


/-I 1 0 0\ 

- 11-1 0 
0 1 - 1-1 
\-l 0 0 0 / 


Since n = 10, we can conclude that each irreducible representation occurs exactly once in the 
isotypic decomposition (5) and therefore the eigenvalues of the adjacency matrix are Ai = 3, 
A 2 = 1, and A 3 = — 2 with multiplicities 1, 5, and 4. 


The irreducible representations were computed with the discrete algebra tool GAP [8], Note that 
the matrices are not orthogonal. Using (7), we compute the projections Pi = , 



/ 

3 

1 

-1 

-1 

1 

1 

-1 

-1 

-1 




1 

3 

1 

-1 

-1 

-1 

1 

-1 

-1 

-1 



-1 

1 

3 

1 

-1 

-1 

-1 

1 

-1 

-1 



-1 

-1 

1 

3 

1 

-1 

-1 

-1 

1 

-1 

1 


1 

-1 

-1 

1 

3 

-1 

-1 

-1 

-1 

1 

6 


1 

-1 

-1 

-1 

-1 

3 

-1 

1 

1 

-1 



-1 

1 

-1 

-1 

-1 

-1 

3 

-1 

1 

1 



-1 

-1 

1 

-1 

-1 

1 

-1 

3 

-1 

1 



-1 

-1 

-1 

1 

-1 

1 

1 

-1 

3 

-1 


V 

-1 

-1 

-1 

-1 

1 

-1 

1 

1 

-1 

3/ 
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P3 = — 
^ 15 


( ^ 

-4 

1 

1 

-4 

-4 

1 

1 

1 

1 \ 

-4 

6 

-4 

1 

1 

1 

-4 

1 

1 

1 

1 

-4 

6 

-4 

1 

1 

1 

-4 

1 

1 

1 

1 

-4 

6 

-4 

1 

1 

1 

-4 

1 

-4 

1 

1 

-4 

6 

1 

1 

1 

1 

-4 

-4 

1 

1 

1 

1 

6 

1 

-4 

-4 

1 

1 

-4 

1 

1 

1 

1 

6 

1 

-4 

-4 

1 

1 

-4 

1 

1 

-4 

1 

6 

1 

-4 

1 

1 

1 

-4 

1 

-4 

-4 

1 

6 

1 

V 1 

1 

1 

1 

-4 

1 

-4 

-4 

1 

6 / 


A 


We have already seen that equivariance properties lead to the existence of multiple eigenvalues 
in a generic way (Proposition 2.8; cf. [1, 4]). For control purposes, it is important to analyze and 
understand this phenomenon in more detail. Here corresponding results by Fassler and Stiefel [6] 
and by Gatermann [9] are particularly useful. They have shown that the block diagonal structure 
induced by the isotypic components can even be refined further by introducing so-called symmetry 
adapted bases. 

Let us explain this result in detail. Let d be an irreducible representation of F. Then we define 
the following projections by 


' ' 7er 


-1 






( 8 ) 


Observe that, if the corresponding unitary representations are used, all the projections defined 
above are orthogonal due to the fact that F C 0(n). 

Theorem 2 . 14 . Let us denote the image of by , /r = 1,... ,n^. Then is an A-invariant 
subspace ofW.Q for every p and, moreover, the restrictions Ay^^ are all isomorphic for p = 1 ,... ,n§. 

An immediate consequence of Theorem 2.14 is the following: Suppose that we make a change 
of coordinates and write A with respect to the symmetry adapted basis. Let d be an absolutely 
irreducible representation of dimension n^. Then there will be identical blocks of dimension 
dim(HA)/n^ on the diagonal of A, and generically these blocks will possess only simple eigenvalues. 

Example 2 . 15 . Let us consider again the system defined in Example 2.11 and compute the sym¬ 
metry adapted basis for the action ofT, see (8). The first two projections remain the same, that is 
Pi = Pi and P 2 = P 2 , and for d^ we compute 


Pf = ^ 


/O 0 

0 



(Id 

0 

-Id 

0 \ 

0 Id 

0 

-Id 

and Po = — 

0 

0 

0 

0 

0 0 

0 

0 

^ 2 

-Id 

0 

Id 

0 

Vo -Id 

0 

Id) 


\ 0 

0 

0 

0 / 
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With this information, we define the symmetry adapted basis 



Wi W2 W3 


and obtain the following block diagonal structure for A by a corresponding coordinate transforma¬ 
tion: 


/22 

-4 

0 

0 

0 

0 

0 

0 \ 

5 

40 

0 

0 

0 

0 

0 

0 

0 

0 

-2 

-16 

0 

0 

0 

0 

0 

0 

1 

20 

0 

0 

0 

0 

0 

0 

0 

0 

10 

-10 

0 

0 

0 

0 

0 

0 

3 

30 

0 

0 

0 

0 

0 

0 

0 

0 

10 

-10 

\o 

0 

0 

0 

0 

0 

3 

30 / 


Example 2.16. For the Petersen graph example, refining these projections with the aid of (8) 
leads to projections P^, P 2 , ■ ■ ■, P|, and P^,..., P^. Here, in order to obtain an orthonormal basis 
we used the corresponding unitary representations. (Note that due to Weyl’s unitary trick 
finite-dimensional representations of finite groups are unitarisable.) The subspaces are spanned 
by the column vectors of the matrix 


(I/Vw 

i/Vw 

1/3 2/y/m 

4/V60 

0 

0 

2/VW 

1/3 

2/Vf^\ 

i/\/Io 

i/\/i8 

1/3 -2/y/m 

i/y/m 

1/2 

0 

-2/y/W 

-1/3 

2/y/m 

i/\/Io 

-l/\/l8 

-1/3 -2/y/m 

l/y/m 

1/2 

l/\/6 

l/\/i8 

1/6 


i/Vw 

- 3 /V 18 

0 0 

0 

0 

-1/V6 

-I/VI 8 

1/3 

2/Vf30 

i/\/To 

-l/\/l8 

1/6 5/\/60 

0 

0 

0 

0 

-1/2 


i/\/Io 

l/\/l8 

-1/6 -l/\/60 

3/\/60 

-1/2 

0 

-2/\/l8 

1/6 


i/\/To 

l/\/l8 

1/3 -2/\/60 

-4/\/60 

0 

-1/V6 

l/\/l8 

1/6 


i/Vw 


-2/3 0 

0 

0 

-l/\/6 

i/\/i8 

-1/3 

2/y/m 


-i/Vw 

1/6 -2>/y/m 

-l/y/m 

-1/2 

1/V6 

i/\/i8 

-1/3 

2/y/m 

\i/Vw 

I/VI 8 

-1/6 3/Vf30 

-4/V60 

0 

1/V6 

- 1 /V 18 

1/3 

2 /V 60 ) 


The matrix T transforms the adjacency matrix into a diagonal matrix whose entries are the eigen¬ 
values ofV. A 
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3 Controllability and Observability in Equivariant Systems 

In this section, we will consider systems of the form 


X = Ax + B u, 
y = C X + D u, 


( 9 ) 


where x G M"' is the state vector, u € the control vector, and y G the output vector. 
Furthermore, A G M”’”, B G MA’P, C G and D G The matrix D is often not required 

and can assumed to be zero. We will introduce only the basic definitions required for our analysis, 
for a more detailed description see for example [10]. 

Controllability and observability of such systems play an important role in modern control system 
theory. Controllability can be defined as follows: If for any given state xq at time to, there exist a 
finite time ti and a control function u{t) that transfers the state x{to) = xq to any state x{ti) = xi, 
then the system is said to be completely controllable. Similarly, a system is defined to be completely 
observable if for any time to the initial state x{to) can be determined by observing the system output 
y{t) over a finite time interval to < t < ti. 


Definition 3.1. Define 

P={B, 

to be the controllability matrix and 


the observability matrix. 


AB, A^B, ••• 

f ^ \ 

CA 

Q= CA^ 

\CA^-^J 


A^-^B) 


( 10 ) 


( 11 ) 


It is well known that a linear time-invariant system of form (9) is completely controllable if and 
only if the matrix P has full rank, i.e. rank(P) = n, and completely observable if and only if Q has 
full rank, i.e. rank((5) = n. These rank conditions are typically used to determine whether a given 
system is controllable and observable. Another way to test for controllability and observability is 
to compute the rank of the matrices 


P = (s/ — A, B) and Q 


si -A 
C 


( 12 ) 


System (9) is completely controllable if and only if rank(P) = n and completely observable if and 
only if rank((5) = n for all s G C. If one of these two equivalent controllability conditions is 
satisfied, then we call the matrix pair [A, B) controllable. Analogously, we call [A, C) observable if 
one of the equivalent observability conditions is satisfied. 

Our considerations will be based on the following well known result (see e.g. [II] or also [4]): 

Lemma 3.2. If {A,B) is controllable and rank(i?) = q, then the geometric multiplicity of each 
eigenvalue of A is at most q. 
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We now assume that A is F-equivariant (see (2)). Let i?!,... ,'i9k be the irreducible representa¬ 
tions of r, denote by their respective dimensions (z = 1 ,..., K), and by P^^ the projection onto 
the corresponding isotypic component (see (7)). For simplicity, we assume that these representa¬ 
tions are absolutely irreducible. The case where F possesses irreducible representations which are 
not absolutely irreducible can be treated in an analogous way. 

By a combination of Lemma 3.2 with Theorem 2.14 we obtain our first result. 

Proposition 3 . 3 . Suppose that {A,B) is controllable with rank(B) = q. Then 


q > A^r, (13) 

where 

Nr = max {m : 7 ^ 0}. (14) 

Proof. Observe that the condition 7 ^ 0 guarantees that the irreducible representation of 
F actually occurs nontrivially in the isotypic decomposition (5). Therefore, by Theorem 2.14 the 
matrix A possesses an eigenvalue of geometric multiplicity at least A^r. Now the result follows with 
Lemma 3.2. □ 


Example 3 . 4 . In our guiding example Nr = 2, and therefore the rank of B has to he at least 
two. For the nonsymmetric case Ci C 2 , however, Nr = 1 and rank(P) > 1 is sufficient. For the 
Petersen graph example Nr = 5, which requires rank(P) > 5 in order to achieve controllability. A 


From now on, we will assume that P^. 7 ^ 0 for i = 1,... ,K. That is, we will restrict our attention 
to the irreducible representations which actually possess nontrivial isotypic components. We will 
now refine the result in Proposition 3.3 and address the question how to use the symmetry adapted 
basis for the identification of appropriate matrices B. That is, we would like to identify a control 
such that both the number of columns of B and the number of nonzero entries of B are as small 
as possible. 

Let us denote by P^., i = 1,..., K , p = 1,... ,ni, the projections onto the symmetry adapted 
basis and by the corresponding subspaces (see Theorem 2.14). Since the projections P^, are 
orthogonal, we know that 

= ( 15 ) 

is an orthogonal sum of the ^-invariant subspaces . We denote by 


lorii nl fon2 nl rW-K 

U)i,...,Di ,02,.. .,02 , . . . , Djr, ■ ■ ■ 1 


(16) 


a corresponding orthonormal basis of M”' where is an orthonormal basis of Vj^, of length di. The 
coordinate transformation induced by this basis is called T. 
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Example 3 . 5 . Considering Example 2.11, we see that 


/ l / 2 \ 


f 0 ^ 


/ l / 2 \ 


f 0 ^ 


f ° ^ 


/ 


0 


1/2 


0 


1/2 


0 



0 

1/2 


0 


- 1/2 


0 


1 /V 2 



0 

0 


1/2 


0 


- 1/2 


0 



l /\/2 

1/2 


0 


1/2 

5 

0 


0 



0 

0 


1/2 


0 


1/2 


0 



0 

1/2 


0 


- 1/2 


0 


- 1 / V 2 



0 

\ 0 y 

\ l / 2 / 


v 0 y 

\- l / 2 / 

✓ ^ 

\ 0 y 

V- 

-1/V2J 


"V* ^ “V* ^ “V 

B\ Bl Bl 


/yV 2 \ 


f ° ^ 

0 


1/V2 

0 


0 

0 


0 

- l /\/2 


0 

0 


-1/V2 

0 


0 

v 0 y 


v 0 y 




yields the orthonormal basis o/M® and thus 


/1/2 

0 

1/2 

0 

0 

0 

1/V2 

0 ^ 

0 

1/2 

0 

1/2 

0 

0 

0 

1/V2 

1/2 

0 

-1/2 

0 

1/V2 

0 

0 

0 

0 

1/2 

0 

-1/2 

0 

1/V2 

0 

0 

1/2 

0 

1/2 

0 

0 

0 

-1/V2 

0 

0 

1/2 

0 

1/2 

0 

0 

0 

-1/V2 

1/2 

0 

-1/2 

0 

-1/V2 

0 

0 

0 

V 0 

1/2 

0 

-1/2 

0 

-1/V2 

0 

0 y 


A 


Proposition 3 . 6 . Suppose that {A, B) is controllable. Then for each irreducible representation Hi 
of r there exist Ui columns 6i,..., of B such that 


P'dfiu ^ 0 for n = 1 ,... ,ni. 
Proof. Consider the block diagonal matrix 

A = P-^ AT = T'^ AT. 


(17) 


(18) 


A necessary condition for (A, B) to be controllable is that for each irreducible representation Hi there 
are Ui columns of B with nontrivial components addressing all the Ui identical blocks corresponding 
to Hi. More precisely, suppose that these n* identical blocks correspond to the index sets I-f C 
{1,... ,n}, p, = 1,... ,nj. Then for each p there is a column b^ of B such that YljeP I 7^ 0- 
Observing that 

B = TB, (19) 

we compute for each such column b^ and 6^ = Tb^ 

= P^Tb^^ = SflP / 0 . ( 20 ) 


Here, denotes the matrix with columns [O -OHf O -O], 


□ 
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Remark 3 . 7 . Consider the generic case, a matrix A with n distinct eigenvalues. For such a 
system, in order to be controllable, the matrix B must also address all eigenspaces. Assume that 
B = there is a j with aj = 0, then 


rank 


\,I-A 


E 


aiVi 


< n 


since Vj G ker{XjI — A) and Vj cannot be generated by the other eigenvectors. Thus {A,B) is not 
controllable. 

With the following proposition we identify particular simple choices for B. 

Proposition 3 . 8 . Recall that di is the length of the orthonormal basis ofVj^, (see ( 16 ) J. Con¬ 
sider the row of the transformation T = and assume tmj / 0 where 


2—1 

j = ^ n^di + (/u - l)di + k 
£=l 


( 21 ) 


for fj. G {1,..., Ui} and a k G {1,..., di}. Then 


pV p 


+ 0 , 


( 22 ) 


where Cm G M"" is the canonical unit vector in M”. 


Proof. The assumption on the index j guarantees that the column Tm of contains a non¬ 
vanishing component corresponding to the invariant subspace in the block diagonal structure 
of ^ = T^AT. Therefore 

P&Tm = P^TT^ = S>fT^ ^ 0 . 


As before, denotes the matrix with columns [0 • • • 0 0 • • • 0], 


□ 


With this proposition, we are now in the position to propose the following algorithm for the 
construction of the matrix B. Let T = (t^s) be the coordinate transformation matrix whose 
column vectors, denoted T^, form a basis. 


1: Set V = 0, define B = B(V) to be the matrix whose columns are the unit vectors 
defined by the indices V. 

2: Choose an absolutely irreducible representation of T with m = Nr, see (14). 

3: Let Ts be the first column of the matrix T belonging to di. 

4: Find the first nonzero entry [Ts]r in the vector Tg, i.e. [Tg],. = Ls A Oj with r ^ V. (23) 
5: V ■(— V U {r} and s ^ s + di. 

6: If \V\ < Ui, go to step 4. (Count only indices added for representation di.) 

7: If rank(P) < n for A and B, choose the next largest representation di and go to 
step 3. 
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In the worst case, the algorithm would add all n unit vectors so that {A,B) is controllable. 
Instead of adding unit vectors directly to the matrix B, one could check whether adding that 
vector would increase the rank of P in order to avoid redundant information in B. Let us now 
apply this algorithm to the guiding example. 

Example 3.9. In this case, Proposition 3.8 and the structure ofT imply that 

P^^es 7 ^ 0 and P^^ei / 0. (24) 

In fact, it can easily be verified that {A, B) is controllable with the choice 


B = 


/I 

0 

0 

0 

0 

0 

0 

Vo 


0 

I 

0 

0 

0 

0 

0 / 


(25) 


Here, Algorithm 23 would start with column 5 ofT, add 6 $ to B, then select column 7, and add ei. 

Since the system is controllable with that choice of B, the algorithm stops. A 

Remark 3.10. Interestingly, for the network with Z 4 symmetry, one vector, e.g. B = (1,0,0,0,0, 0, 0, 0)^, 
suffices to render the system controllable. That is, the system is controllable from any node. 

Example 3.11. Similarly, using Algorithm 23, we obtain for the Petersen graph example and the 
representation 

^ -^ 12®2 / 0 , 0 , 0 , Pj^eg 7 ^ 0 . 

The corresponding matrix 

/I 0 0 0 0 \ 

0 10 0 0 
0 0 10 0 
0 0 0 0 0 

0 0 0 0 0 

0 0 0 1 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 1 

Vo 0 0 0 0 / 

also addresses the other blocks corresponding to Hi and H 3 . Thus, {A,B) is controllable. Not all 
configurations with N-p = 5 input vectors are controllable. Figure 3 shows different configurations 
which are controllable (controlled vertices are marked green) and noncontrollable (controlled vertices 
are marked red). A 
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Figure 3 : Controllable networks (in green) and noncontrollable networks (in red) with five inputs. 


So far, we focused only on the controllability of a given system, but due to the duality of 
observability and controllability - (^4, C) is observable if and only if , C'^) is controllable -, all 
the results described in this section can also be used to determine observability. 


4 Conclusion 


In this paper, we analyzed the relationship between symmetries in complex networked systems and 
the controllability and observability of such systems using results from representation theory and 
illustrated that symmetries, depending on the symmetry type, might decrease the controllability. 
Furthermore, we derived necessary conditions on the number of inputs to guarantee controllability 
and presented an algorithm which computes a sparse input matrix B based on projections onto the 
isotypic components. We demonstrated the algorithm using two guiding examples, a network with 
symmetry and the Petersen graph, whose automorphism group is isomorphic to S^. 

Future work would be to analyze how the results presented within this paper can be extended 
to time-varying systems of the form 


X = A{t) X + B{t) u, 
y = C{t) X + D{t) u. 


( 26 ) 


For such a system, controllability and observability can be defined in terms of the so-called Control¬ 
lability Gramian and Observability Gramian. One possible extension would be to analyze structural 
controllability and observability of these networks using the nonzero patterns of the matrices. The 
time-varying network structure could result in a temporarily uncontrollable system if a matrix 
A(fo) without symmetries at time to is changed into a highly symmetric matrix A{ti) at time ti. 
Our results show that for such a system the number of required control inputs might increase 
significantly. Consider for instance the guiding example: If Ci ^ C 2 , the system is controllable 
from one node, for Ci = C 2 , on the other hand, two control inputs are required. Equivalently, 
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adding links to or removing links from a network might cause new symmetries or break existing 
symmetries and, as a result, change the number of required control inputs. 
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